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I extend the usual linear-theory formula for large-scale clustering in redshift-space to include 
gravitational redshift. The extra contribution to the standard galaxy power spectrum is suppressed 
by k~ 2 , where k c — ck/aH (k is the wavevector, a the expansion factor, and H = a/a), and is thus 
effectively limited to the few largest-scale modes and very difficult to detect; however, a correlation, 
tx fc" 1 , is generated between the real and imaginary parts of the Fourier space density fields of 
0^ ■ two different types of galaxy, which would otherwise be zero, i.e., the cross-power spectrum has an 

imaginary part: P ab (k, /J,)/P(k) = (b a + ffi 2 ) (b b + / M 2 ) - i§fi m £ (b a - h) + 0(k~ 2 ), where P(k) 
is the real-space mass-density power spectrum, bi are the galaxy biases, \x is the cosine of the angle 
between the wavevector and line of sight, and / = dlnD/dlna (D is the linear growth factor). 
The total signal-to-noise of measurements of this effect is not dominated by the largest scales - it 
converges at k ~ 0.05 ftMpc . This gravitational redshift result is pedagogically interesting, but 
naive in that it is gauge dependent and there are other effects of similar form and size, related to the 
0^ ' transformation between observable and proper coordinates. I include these effects, which add other 

contributions to the coefficient of n/k c , and add a fi 3 /k c term, but don't qualitatively change the 
picture. The leading source of noise in the measurement is galaxy shot-noise, not sample variance, 
so developments that allow higher S/N surveys can make this measurement powerful, although it 
U ' would otherwise l>o only marginally detectable in a -IDOI-seale survey. 
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Redshift-space distortions by peculiar velocities were first computed in [I], leading to the well-known formula for 
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I. INTRODUCTION 



the galaxy power spectrum, P g (k,fi) = (b + //i 2 )" P(/c) where 6 is the bias, / = dlnD/dlna, D is the linear growth 
factor, a the expansion factor, fi is the cosine of the angle between the wavevector k and the line of sight, and P(k) 
is the real space mass-density power spectrum. Measurements of these distortions can provide powerful constraints 
on cosmology, isolating the bias-independent quantity f 2 P(k), including it's redshift evolution [2, 0, 0, HI, 111 • It has 
^v^i ' always been known that the formula of [l[ is approximate in a variety of ways, which must be improved as precision of 
^vq , measurements improves. The most important corrections are for non-linearity, which are important on relatively small 
scales where we have the most statistical power @, H, H, E3 EH EH EH • The form of redshift-space distortions can 
be modified even on linear scales by non-linear transformations of the alre ady distorted field [3, EE Ell or possibly 
selection effects [l7j]- Finally, on large scales there are relativistic effects O Eil, which will be the subject of this 
Q\ ; paper. 

We normally assume that the Universe is intrinsically homogeneous and isotropic. Our line of sight breaks this 
symmetry when we perform a redshift survey, however, standard redshift-space distortions due to peculiar velocities do 
not break reflection symmetry along the line of sight, i.e., the survey would look the same if viewed from the opposite 
' direction (ignoring effects on the photons as they travel to us, and anything else not included in the standard result 
of [H) ■ Gravitational redshifts do break this symmetry, e.g., the redshift of a photon coming out of a potential well 
is independent of viewing angle, so the apparent displacement is always away from the observer, i.e., it changes sign 
if the viewer moves to the opposite side. Note, however, that when measuring a standard auto-correlation function 

£ m| = X|| — ^ = ^<5(x||)(5(x'|)^ (where x\\ is the radial coordinate and I have suppressed the transverse coordinates 

because they are irrelevant to this discussion), there is nothing to distinguish £(r||) from £(— rii), i.e., if the galaxies 
are identical, there is no way to detect an exchange of their positions. If we consider the cross-correlation of two 

types of galaxies, a and 6, the situation changes - there can in principle be a difference between ^<5 a (a;|| )Sf,(x'^ )\ and 

(^5b(x\\)5 a (x^) \ , i.e., equivalently, between £ a b(f||) and £ a &( — r \\)' The power spectrum, the Fourier transform of the 
correlation function, is normally real because the correlation function is even in . If the correlation function has 
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a component that is odd in rn, the power spectrum has an imaginary part. Such an imaginary part, generated by 
gravitational redshift and other redshift-space distortions not included in will be the subject of this paper. 

The plan of the paper is as follows: In $H]l compute the effect of gravitational redshift alone. This calculation is 
intended to be pedagogical, presenting the basic form of the new effect in a calculation that is easy to understand - it 
does not contain all of the terms of similar magnitude and is not gauge invariant. In ^IIII I add the additional terms 
recently computed by [l9j |. and estimate how precisely it will be possible to measure this effect in the future. Finally, 
in ^IVU briefly summarize and discuss the results. 



II. NAIVE GRAVITATIONAL REDSHIFT EFFECT 



I first compute the redshift-space distortion due to gravitational redshift in the way that standard velocity-induced 
distortions have been computed It turns out that other terms are similarly important, but this calculation contains 
all of the qualitatively new k dependence (where k is the wavevector). I am going to compute what for the CMB 
would be called the primary fluctuations, dependent only on the potential at the location of the galaxy, not on an 
integral along the photon's path to us. 

In the small-separation limit, the radial coordinate measured by a redshift survey is 



cAA 



HaAxu + Avu 



c 



(1) 



where the first term is the usual Hubble expansion of the Universe (H — a /a), the second is the Doppler shift due 
to radial peculiar velocity vn, and the last is gravitational redshift due to potential ip. The gravitational redshift 
is usually ignored, for what we will see has been a good reason, for relatively small-volume existing surveys. The 
apparent comoving distance is 



As|| = Aa 



(Hay 1 Av\\ - (cHay 1 Atp 



(2) 



I assume that the conversion from transverse angular separation to Hubble velocity separation is known, so the 
components of x and s transverse to the line of sight are equivalent. Including the uncertainty in this angle- velocity 
conversion leads to the Alcock-Paczynski test. In addition to assuming that we know the functions H(z) and D^(z), 
I also assume that we know at what redshift to evaluate them - of course, for a real survey we do not, we only know 
the observed redshift of the galaxies, not the expansion factor in the background Universe. This, and similar issues, 
will lead to the additional terms, computed by [Tj|, which I add in the next section. Density in redshift space, p s is 
related to density in real space, p, by p s d 3 s = pel 3 *.. Using 
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and making the distant observer approximation [20(, gives 

(1 + S 8 ) 



l + (Hay^-(cHay^ 

OX\\ " x \\. 
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= (1 + 6(5) 

dip 



S s =bS- (Hay 1 ^1 + (cHay 1 Q . 

where I am working to linear order in the perturbations. In Fourier space this is 

S s k = b <5k + i (Hay 1 fc||f||k - i (cHay 1 k\\i(fk ■ 



The velocity is given in linear theory by 
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(4) 



(5) 



(6) 
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where p — ku/k and k c — c k/a H is the wavevector in the observable Hubble velocity units over the speed of light 
(v is also normalized by the speed of light). The potential is 



fl m (a) k c 2 (5 k 



(8) 
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Then, 



(9) 



Q m — fl m (a) here, and in general such quantities should be assumed to be time dependent, rather than the z — 
value, unless otherwise indicated. The power spectrum of (5 s k is 



(S sk S sk ,) = (2irf 5 U (k + k') 



(6 + //i 2 



2 m kr 



P(k) 



(10) 



where -P(fc) is the real space mass density power spectrum. The fact that the new term is proportional to (ck / aH)~ 2 
means that it is only significant on very large scales, if at all. It is useful to plug in some numbers to see just 
how bad this is. A 100 cubic Gpc/h JDEM-scale survey [2l| (~ 2/3 of the sky over 1 < z < 2) has a minimum 
k ~ 2n/V 1 ^ 3 — 0.0014 h Mpc -1 or k c ~ 4.5, i.e., one should not be tempted to think that k c ~ 1 is reachable 
(every bit of volume out to z ~ 5 would only reach k c ~ 2). This means that, for the 100 cubic Gpc/h survey, 

for the best case fi = 1, the quantity ^f^m-^-) = 0.08 for the largest (fully sampled) modes in the survey. The 

standard power, which sets the noise level, has (b + ffi 2 ) 2 = 3.7, for a modest 6=1, and f(z = 1.5) = 0.92, i.e., 45 
times larger than the gravitational redshift effect, which clearly cannot be detected in the straightforward auto-power 
spectrum of galaxies (the k~ 2 decline in the signal means that one can not sum many modes to overcome the small 
size of the effect). On the bright side, if one is interested in detecting non-Gaussianity in the LSS power spectrum 
[H, [H, [H, HE H3, [H, HI, [33, HH, SIM H, [H , the smallness of this kind of effect is a good thing (there will be 
other, isotropic, terms of similar order [19j |. which would look like the effect of the local model of non-Gaussianity). 

We would really like to find a cross-term between the gravitational redshift term in Eq. © and the standard terms, 
which would be proportional to k^ 1 instead of k~ 2 , and we would also like to somehow avoid competing with the 
standard power as background noise. It turns out that if we cross-correlate two different types of fields with different 
bias we can satisfy both of these desires. 



(6^,) = (2nf5» (k + k') 
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(11) 



Note that one can isolate the imaginary term operationally by computing 

Pf(k) = (Im [6^6%]) = (6%S% 
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where S l s ^ and S l J k are the real and imaginary parts of the galaxy density field. This is directly computable from 
observations, without knowing anything about the underlying theory. We generally would not compute it because it 
is assumed to be zero. 

It may not be completely obvious how to compute the expected error on a measurement of Pj b (k), and, for the 
next section, it will be useful to go through the whole calculation for a galaxy density field of the following general 
form: 



(13) 



where 5^ and 5^ are the real and imaginary parts of the mass density field and i?^ and are the real and imaginary 
parts of the bias- type coefficients, for galaxy type i (i.e., / would traditionally be zero, while in the calculation of this 
section Ik = |fi m #-). Then 



(im [«j° k <&]) = (4X - Kil) (SL 2 + s5 2 ) = (4X - Riit) P(k) 



(14) 



A key fact about this calculation is that, if P = 0, Im [<5" k ^sk] = mode-by-mode, i.e., before taking any expectation 
value. This means that the sample variance will also go to zero if I — > 0. Note that I have been ignoring shot-noise. 
If I continue to ignore it, I find the variance 



Im [5 



(lm[5 s V s 6 k ]» 2 ) = (J^ 



: P(fc) 2 = (lm[^]) 2 



(15) 



i.e., the error on (im [<5" k <5*£] ) follows the usual rule for a power spectrum, with S/N = 1 per mode. Remarkably, 
this shows that the fractional errors on the imaginary part of the power spectrum in the cosmic variance limit are 
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just as small as the fractional errors on the standard power spectrum, in spite of the fact that the real part is much 
larger and one might have expected it to provide background noise. This feature relies on the fact that the fields for 
the two types of galaxy trace the same underlying fluctuations. Once noise is added, the situation is not quite as 
rosy, although still quite good. I assume the noise is the standard Poisson sampling noise, uncorrelated between the 
two types of galaxy. I write the measured density as <5* k = <5* k + e k . The noise makes no mean contribution to the 
imaginary part of the cross-power spectrum, but it does add variance in the measurement: 



Im 



°sk°sk 



P 7 ab (k)M = P? 6 (k) 2 + ^N a P b (k) + ^N b P a (k) + ^N a N b , (16) 



where N l is the noise power for galaxy type i, and P l (k) is the auto-power spectrum for galaxy type i (always real). 
In the relevant regime, the errors here will be dominated by the terms containing P l cx R 2 . Remember, however, 
that the imaginary signal itself is proportional to RI, so, at an order of magnitude level (assuming the biases do not 
nearly cancel), the i?'s cancel in the S/N ratio, leaving (S/N) 2 oc I 2 P/N, i.e., the detectability of the signal depends 
on the size of the gravitational redshift effect compared to the noise, not compared to the real part of the signal. As 
the noise is reduced, the measurement improves, until one reaches the limit where S/N = 1 per mode (in practice, 
this limit will be very difficult to reach). More concretely, 



S\ 2 ^2(l-Rl-R-I b k ) 2 P(k) 



Nj k R b2 N a + RfN b 



(17) 



The signal-to-noise of a power spectrum measurement aimed at a simple detection of some effect in a survey with 
volume V is roughly 



, , k 2 dk / dul — (k,u) . (18) 



If (S/N) 2 (k) (x k 1 with 7 < —3, the integral converges quickly at the low k end, while if 7 > —3 the total S/N 
continues to increase with increasing A; max . For S/N(k) rolling with k, the integral will be dominated by the scale 
where 7 ~ —3. In the case at hand, j(k) = n c g(k) — 2, where n c s(k) = dlnP/dhik, i.e., the integral is dominated by 
k ~ 0.05 h Mpc" 1 , where n c g ~ —1. I will leave a numerical evaluation of the detectability of the signal for the next 
section, where I include several other effects which change the signal amplitude in detail, although not the order of 
magnitude or form. 



III. THE REDSHIFT-SPACE POWER SPECTRUM TO O^" 1 ) 

Recently, presented a much more rigorous and complete calculation of the observable redshift-space galaxy 
density. Here I will use their result to compute the redshift-space power spectrum, which amounts to changing the 
coefficients in I k above. When we start to include relativistic terms, it becomes natural to look at the calculation 
as an expansion in l/k c . I will compute the first, 0(k~ 1 ), correction to the standard LSS picture, dropping terms 
0(k~ 2 ), whi ch as I discussed above, are very hard to detect (except in projection against the very smooth CMB 
background HEHEHl). 

Before considering galaxy clustering, it is useful to check if there are any k~ x corrections to the standard calculations 
of the perturbations in mass-density, 5, velocity, v l , and potential, ip, fields. I will consider the evolution equations 
appropriate for dark matter perturbations only, because we are interested in late times and large scales where baryons 
and dark matter are equivalent, and radiation is negligible (I assume a cosmological constant, so no dark energy 
density fluctuations). I will present the results in terms of quantities calculated in the conformal Newtonian gauge, 
where the metric is: 

ds 2 = g^dx^dx" = a 2 (r) [- (1 + 2%j)) dr 2 + (1 + 20) dx 2 ] . (19) 

-0 can be identified with the Newtonian potential, and <f> — —tp. The final galaxy clustering results will be gauge 
invariant, because they are calculated entirely in terms of observable combinations of the gauge-dependent quantities. 
The evolution equations are [39} : 



5 + 30 + e = , 



(20) 
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where & = diV % and the dot is a derivative with respect to conformal time (note that this is the linearized Newtonian 
continuity equation plus an extra <f> term), 

6 + HO - k 2 iP = (21) 
where 7i = Ha (just the standard linear Euler equation), and 



4irGa 2 p„ 



m 

k 2 



(22) 



(the Poisson equation, plus the extra 9 term). We see immediately that all corrections are 0(k~ 2 ), i.e., irrelevant 
to this calculation. We could have guessed this from symmetry - all of the effects at 0(k~ 1 ) will be related to the 
observing process, where reflection symmetry is broken along the line of sight . 

Most of the work needed to compute the galaxy power spectrum has been done by [li^ |. Their equation Eq. (36), 
which I write here in terms of conformal Newtonian gauge sub-components, gives the gauge invariant fluctuations in 
galaxy density, that one can observe using only the measurable frequency of photons and angular position on the sky: 

„ .. , dS z 1 + 2: „ „ 1 + z dH 5 r , 

Sobs = b(6 m ~ 3(5, + V + 20 + uii - 1 + z)—^ - 2——S z - 6 Z - b P 8v L -2n+ -j^—^S, + 2^ , 23 

" Oz Hr H dz r 

where the various elements of this equation, discussed below, are derived in [l^ | . I have not verified that this equation 
is correct. Note that it is certainly not perfectly complete, e.g., evolution in the number density of galaxies would 
produce at least one more term (this would not induce any interesting new behavior, so I continue to neglect it). 
The redshift z is the observed one, and background quantities like H are evaluated at this observed z. I am going 
to assume that the background quantities are known functions of z. Relaxing this assumption leads to the usual 
Alcock-Paczynski-type effects |4fJ, which, it should be understood, are different in nature from the ones that lead 
to Eq. (|23| . The terms in Eq. (|23|) are related to the uncertainty in the z at which to evaluate the background 
functions (among other things), rather than uncertainty in the functions themselves. One can perfectly well imagine 
studying LSS in a Universe where the background evolution was known perfectly (e.g., for the purpose of studying 
the primordial perturbations), and in that case Eq. (|2"3"]) would be the one to use. 

I now pare Eq. ([231) down to 0{k~ l ) terms, in addition to dropping the k (lensing) term for simplicity. Lensing 
can be important [4l| , but the fact that it is sensitive to perturbations integrated along the line of sight means that it 
produces a power spectrum with a much different form than the terms that I consider, outside the main point of this 
paper. Everything I consider here would be called primary fluctuations in CMB language, i.e., not depending on an 
integral over perturbations along the line of sight. Eq. (16) of [l9| shows that 5r/r ~ 0(k~ 2 ) so it can be dropped. 
Note that r is the comoving distance to the galaxy. The redshift perturbation 8 Z can be inferred from Eq. (11) of 
[T^ |. and the subsequent paragraph: 



ip- dr (i>~ <p) (24) 



or 

S z = v\\ + 0(k~ 2 ) (25) 

and 

f^-§W c - 2 ) (26) 
oz oz oz 

(Recall that v oc k~ 1 5, ip oc k~ 2 6, and a derivative brings in an extra factor of k.) 
Eq. (27) of [H gives 

1 + *, . ~„_ 2 



Hr 



5 Z + 0(k-') (27) 



where I have dropped the lensing component. Note that the luminosity bias term should be viewed as a place-holder 
for a variety of possible similar terms that can appear, depending on the type of observable, i.e., it assumes our 
observable is a simple observed-magnitude-limited galaxy density, while in reality we may be using other things like 
halo- mass- weighted density [42[ or 21cm intensity mapping [43| or the Lya forest [44114^. One might be tempted to 
make a "distant observer" approximation and drop the term oc c/TLr (note that I usually use c = 1), however this 
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is not a good approximation, e.g., at z = 1.5, c/TLr ~ 1. In the language of this paper, the usual justification for 
dropping this kind of term is actually the k~ x expansion, only breaking down if 1/r becomes large enough to overcome 
this. 

Using (1 + z)d z =TL~ 1 d r , we are left with 
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lensing 



(28) 



We see that the first two terms are the usual redshift-space density field, the third is the gravitational redshift term 
computed in CTfl and then there are a variety of terms proportional to vn. 

Moving to Fourier space, we see that the v\\ term has the same k dependence as the gravitational redshift term, 

i.e., vjj. = —if^-Sic, while —ik\\ipk/H = i|O m ^-5k, the difference is only in the coefficients. Finally, 
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(29) 



Note that the Fourier transform in the radial direction involves a bit of a slight of hand, as the statistics of the field 
are not stationary (distance-independent) in that direction. This is not a new problem, as similar redshift evolution 
has always been present in this kind of calculation within the growth factor (hidden in 5(z)) and f(z) and b(z), but 
it is highlighted by the explicit presence of the comoving distance r(z) itself (which is also the coordinate in the 
Fourier transform (it may be more comforting to think of Hr as the Alcock-Paczynski factor)). Of course, the Fourier 
transform in this type of calculation is also somewhat formal because the surveys are not infinite or periodic. Dealing 
with these issues will require very careful data analysis, presumably done at least partially in configuration space. In 
the end, although they will make one's life more complicated, it is likely that background uncertainty and evolution 
issues will lead to new opportunities to derive information from these surveys, as exemplified by the Alcock-Paczynski 

test m. 

I now use Eq. (I14|) to compute the imaginary part of the power spectrum for two types of galaxy, 
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2 m J \ Hr 
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(30) 



One thing to note about this equation is that 



2 r 



+ f 



dlnH 
line 



in an Einstein-de Sitter Universe, and will therefore 



be fairly small in the real Universe, except at low z (where there is not a lot of volume). 1 — ^- is not necessarily 
quite as generally small, but is nearly zero at z ~ 1.5 where a JDEM-like survey would be targeted. Beyond that, 
it is useful to plug in some numbers in order to do a concrete calculation of detectability. I will assume b a = 2, 
bb = 1 (this level of difference is quite reasonable, e.g., the BigBOSS proposal [46[ contains two different samples with 
rou ghly these biases, although not observed on quite the scale assumed here), p a = 1, Pb = 0.5 (roughly reasonable 
[HS3, and not very important), and tt m . = 0.28 This gives (ftr)- 1 = 1.08 at z = 1.5, n m (z = 1.5) = 0.86, 
f(z = 1.5) = 0.92, and dlnH/dlna = -1.29. This makes Pf b {k,p) = {p/k c ) (-2.1 + 2.7V 2 ) P(k). Finally, I assume 
a 100 cubic Gpc/h survey, and noise given by n l P l {k = 0.2 h Mpc -1 ) = 1, where n l = \/N l is the number density 
of galaxies of type i. This noise level is often used as a target for BAO surveys, although nP ~ 3 is really required 
to squeeze out most of the BAO information. Integrating Eq. (fTH)) . I find S/N — 3, when the integration is taken to 
k = 0.2 /iMpc -1 (S/N = 2.7 if the limit is 0.1 /iMpc^ 1 ). To scale this to other surveys, S/N cx V 1/2 and S/N oc n 1 / 2 , 
i.e., we would reach a 5-rr detection for nP — 3. 



IV. CONCLUSIONS 



The main result of this paper is Eq. (f3"0"]) . the imaginary part of the cross-power spectrum between two different 
types of galaxies, which is generated by gravitational redshift and other effects of similar form in the redshift-space 
density field. The signal is suppressed by ^ relative to the standard power spectrum, but the increasing number of 
modes at higher k makes the typical scale contributing to a detection relatively small, i.e., k ~ 0.05 /iMpc" 1 . This 
effect should be detectable in JDEM-scale [2l[ surveys, although it does require subsets of galaxies with substantially 
different bias, both with good signal power-to-noise power ratios, so a bare-bones BAO experiment is not necessarily 
optimal - a 100 cubic Gpc/h survey that can achieve nP(k — 0.2 /iMpc -1 ) = 1 for two samples of galaxies with a 
separation in biases of ~ 1 would make a ~ 3 — a detection. The measurement errors will be dominated by shot-noise, 
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not sample variance, so there is a lot of room for improvement. 21cm intensity mapping [43| could provide a perfect 
low bias, high S/N field complementing a high bias B AO-oriented galaxy survey. Creative halo weighting schemes 
[42l ] promise to change the way we think about shot-noise, possibly reducing the noise power in planned surveys by 
an order of magnitude below n^ 1 . 

Ta king the imaginary part of the power spectrum naturally has the effect of beating cosmic variance in the sense 
of 0, l49j ] so it isn't clear that anything more can be gained in that direction. This result continues the recent trend 
toward viewing the fact that different kinds of galaxies have different biases, i.e., trace the density field in different 
ways, as an opportunity to tease out different physical effects, rather than simply as a nuisance [2j, |42j, |49( . While the 
short-term planning of surveys will probably continue to follow the "measure the BAO feature by looking for high 
bias and volume with nP ~ 1" strategy, it is clear that the long-term future of LSS studies is much richer than that. 
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